Nature of the antiferromagnetic to valence-bond-solid quantum phase transition 
in a 2D XY-model with four-site interactions 
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We report large-scale quantum Monte Carlo calculations at the T = antiferromagnetic to 
valence-bond-solid (VBS) transition of a two-dimensional S — 1/2 XY model with four-spin in- 
teractions. Finite-size scaling suggests a discontinuous spin stiffness, but a continuous VBS order- 
parameter. We propose that this is a continuous VBS transition without critical spin fluctuations. 
We also argue that the system is close to a point, in an extended parameter space, where both the 
magnetic and VBS fluctuations are critical — possibly a deconfined quantum-critical point. 
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Senthil et al. recently proposed a class of generic decon- 
fined quantum- critical points describing phase transitions 
between 0{2) or 0(3) antiferromagnetic (AF) and four- 
fold degenerate valence-bond-solid (VBS) ground states 
in two dimensions (2D) |lj. These critical points fall out- 
side the standard Landau-Ginzburg- Wilson framework, 
where order-order transitions are generically first-order 
(except at fine-tuned multi-critical points). A continuum 
field theory of spinous interacting with a U{1) gauge field 
was proposed in which the critical point is not explicitly 
described in terms of order parameters. Instead, AF or 
VBS order is a consequence of confinement of spinous. 
This remarkable, but so far untested, theory calls for nu- 
merical studies of lattice models that could potentially 
exhibit deconfined quantum-criticality. It was suggested 
[ij that such a transition may already have been ob- 
served in a quantum Monte Carlo (QMC) study |3l of an 
XY-model including four-site interactions. Such a model 
should also from a theoretical perspective be among the 
most natural candidates for the proposed physics [ij, |3| . 

The J-K model studied in 2] is defined by 
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of the J-K model in the vicinity of the AF-VBS tran- 
sition. We have performed ground-state (T — *■ con- 
verged) finite-size scaling for L x L lattices with L up to 
K, 100. Although the spin stiffness shows signs of a dis- 
continuity, we argue that the transition is not necessarily 
first-order. We propose a scenario in which VBS order 
can emerge continuously in a state where AF order is 
not suppressed by diverging VBS fluctuations. The spin 
correlations change continuously from long-ranged to ex- 
ponentially decaying at the transition, while the stiffness 
is discontinuous because the VBS formation is associ- 
ated with the opening of a spin gap. While our results 
show that the theory of deconfined quantum-criticality 
does not apply, we argue that such a transition may be 
located nearby in an extended parameter space. 

The spin stiffness ps is calculated with the the SSE 
method in the standard way [a, |3 in terms of winding 
number fluctuations. At a quantum phase transition with 
dynamic exponent z, it should scale with size as [Sj 
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We also compute the squared VBS order parameter {rn^p) 
and the associated susceptibility xp\ 



and Pijfcz are, respectively, operators acting on ("^p) — TF7 7 ,{PaPb){—'^Y'' 



where Bi 

nearest-neighbor sites and four sites on the corners of a 

plaquette on a 2D square lattice; 
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This S — 1/2 model is equivalent to a half-filled lattice of 
hard-core bosons. It was found to undergo a continuous 
AF-VBS transition at K/J « 7.9, with no intervening 
coexistence region. However, recent simulations of re- 
lated bosonic current models have instead shown weakly 
first-order transitions 1^ . This motivates us to revisit the 
J-K model and analyze the transition in greater detail. 

Using a stochastic series expansion (SSE) QMC algo- 
rithm |5l lal, we have carried out extensive simulations 



Here Pa — Pijki and the sums are over all plaquettes. 
The expected quantum-critical scaling is [SJ 
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where 77 is the correlation function exponent. 

Analyzing the stiffness according to Eq. Q), we ob- 
tain z « 0.4. In Fig. n^a) we show how L^'^ps versus 
K/J graphed for different lattice sizes produces a cross- 
ing point at K/J « 7.915 (with z — \ the crossing points 
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FIG. 1: (Color online) (a) Dependence of the size-scaled stifi- 
ness on the coupling K/J of the J-K model, using z = 0.4. (b) 
The size-scaled stiffness of the anisotropic Heisenberg model 
versus the Ising anisotropy A, using z = 0.3. 



move significantly with L). As we will show below, this 
unusual value of the dynamic exponent is not consistent 
with the r > behavior of the spin susceptibility. The 
discrepancy could be interpreted as a quantum-critical 
point violating hyperscaling, but a more plausible sce- 
nario is that the observed scaling is an artifact of the 
limited range of system sizes available. We therefore ex- 
plore an alternative scenario. 

Consider the anisotropic 2D Heisenberg model; 
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Its xy spin stiffness is non-zero at T = for A < 1 
and vanishes in the thermodynamic limit for A > 1. 
This is not due to a phase transition, but a consequence 
of the order parameter flipping from the xy-plane to 
the 2:-axis. Exactly at the isotropic point the stiffness 
should approach a constant value; 2/3 of the stiffness of 
the symmetry-broken state (reflecting rotational averag- 
ing). However, numerical results do not show a point at 
which ps becomes obviously size independent. Instead, 
as shown in Fig.^b), we find curve crossings for L'^'^ps 
very close to A = 1, reminiscent of the results for the 
J-K model in Fig. QJa). This shows that a discontinu- 
ity in the thermodynamic limit can easily be mistaken 
for a continuous transition with an anomalously small 
z — going to very large lattices we would eventually find 
z = 0. We therefore believe that the stiffness of the J-K 
model is also discontinuous. This does not have to imply 
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FIG. 2: (Color online) Temperature dependence of the spin 
susceptibility for values of K/J in the vicinity of the T = 
AF-VBS transition. Error bars are not shown but are at most 
the same size as the symbols. 



a point with enhanced symmetry. A similar behavior can 
be expected also if Ps is discontinuous for other reason 
(except in the case of a real level crossing). 

The dynamic exponent can also be extracted from the 
temperature dependence of the uniform susceptibility; 
Xu = JT-^{iT,tS!)^)/N. For a 2D quantum-critical 
system it should scale as |2( 
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where 6 is a constant related to the spin-wave velocity 
and a = at the critical coupling. Away from the critical 
coupling a ^ and there is low-T cross-over to a different 
form. The earlier results |2| were consistent with a ^ = 1 
quantum-critical point; Xu oc T at Kc/J « 7.91. We now 
have higher precision at lower temperatures. Results for 
L = 256 (infinite-size converged) are shown in Fig. |2 
We observe a linear behavior over an extended range of 
temperatures, but a line fitted to the data has a small 
but clearly non-zero intercept, suggesting AF order up 
to Kc/J — 7.915 ± 0.005 [note the consistency with the 
stiffness crossing point in Fig. ^a)]. Thus the AF order- 
parameter is not critical at the transition, in accord with 
a stiffness jump. However, the linearity of Xu{T) sug- 
gests that the system is close to a quantum- critical point 
with z — 1 in an extended parameter space. The devia- 
tions from linearity seen in Fig. |^ at low temperature are 
consistent with gapped spin fluctuations (exponentially 
decaying x«) for K > Kc and an expected |l(| cross-over 
to T-independence for K < Kc- 

Considering the discontinuities we found above, one 
might conclude that the transition should be first-order. 
We do not find any signs of this in the VBS correla- 
tions, however, and will therefore consider the possibility 
of critical VBS fluctuations. Fig.Olshows log-log plots of 
(TTip) and xp versus L in the vicinity of the transition. 
The behaviors are typical of quantum-critical scaling. 
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FIG. 3: (Color online) Finite-size scaling of the squared VBS 
order parameter (a), and susceptibility (b). In (a) the line 
has slope —1, and the ones in (b) have slopes 1/2 and 1. 



From the slope of (m'j,) we get z + 77 « 0.98 ± 0.04 (tak- 
ing into account our estimated accuracy of Kc) using l(7|l. 
The uncertainty in 77 extracted from (jH)) is much higher, 
however, due to the significant changes in xp close to 
Kc- We can only give a rough estimate, rj g —(0.5, 1.0), 
which gives z £ (1.5,2.0). Since z is normally integer, 
this would suggest z — 2, but clearly further work is 
needed to confirm this. A negative rj is unusual, but the 
combination z -I- ry « 1, corresponding to '--^ 1/r VBS cor- 
relations, is not unusual. We note that unusually large 
lattices, L « 40, are required before the (likely) asymp- 
totic behavior of (ttt-p) and xp commences. This could 
be explained as a cross-over due to another critical point, 
located nearby in an extended parameter space, which we 
have already argued for above. For a first-order transi- 
tion, one might expect such a cross-over for K < Kc to be 
followed by a rather sharp change to an 1/L^ behavior for 
both (wp) and xp (reflecting short-range correlations), 
because the cross-over should then correspond the size of 
a typical VBS droplet. Instead, (wp) for both K < Kc 
and K > Kc approaches an « 1/L scaling, and the drop 
in Xp occurs only for much larger L, well inside the AF 
phase. Although we can still not completely exclude a 
first-order transition, the intriguing possibility of critical 
VBS fluctuations but discontinuities in the spin sector at 
the same point deserves further analysis 

A discontinuous ps can be accounted for if the VBS 
would be associated with a continuously opening spin 
gap, as illustrated in Fig. 2Ia) — the stiffness must van- 
ish once an infinitesimal gap opens. In this scenario the 



spin correlations would still change continuously, in a 
non-critical manner, from long-ranged to exponentially 
decaying. To argue for the possibility of such an exotic 
transition, we first discuss the conditions for symmetry 
breaking in the AF and VBS phases in terms of quantum 
levels for finite system size A^. We then assume a continu- 
ous onset of VBS order and investigate the consequences 
of this on the relevant quantum states. 

AF order is associated with a "tower" of quantum-rotor 
states (global spin-rotation excitations) that become de- 
generate with the ground state as 1/A^ when N ^ 00. In 
a spin-isotropic system the ground state has spin S — 
and the rotor states S = 1,2,... ^. The XY model 
does not conserve total spin, but we will use the notation 
S — 0,1,... also for its rotor states. For simplicity we 
restrict the discussion to the S'^ = sector. The ground 
state has momentum fc = and the A: > spin wave 
excitations all have their own rotor towers. 

A columnar VBS state is a quadruplet corresponding 
to the ^4 symmetry of a columnar dimer pattern. We 
refer to the lattice-symmetry related quantum numbers 
with a single label, p = 0, 1, 2, 3, and call the states with 
p — 1,2,3 VBS states. Like the ground state {p = 0), 
the VBS states have S = 0, and, due to the discrete VBS 
order parameter, the level spacing within the quadruplet 
vanishes exponentially with TV. 

We now explore an AF-VBS transition involving only 
the states discussed above, i.e., those states are assumed 
to stay at low energy and no other states descend from 
higher energies. We will follow the evolution of the levels 
\S)p (5 = 0, 1, . . .; p = 0, 1, 2, 3) as we vary K. Although 
the character of the states changes as they evolve, we 
continue to call them rotor and VBS states also in their 
"wrong" phase. We assume that the ground state |0)o 
evolves continuously and is not crossed by any other level. 

For S = Kc — K > the VBS correlation length grows 
as S'" , corresponding to a gap A ~ S^'^ between |0)o and 
|0)i. Thus, for all K < Kc there is some N above which 
there are rotor states |1, 2, . . .)o below the first VBS state 
|0)i. As iV ^ 00 a tower builds up and AF order can 
form. There is no apparent reason why the AF order has 
to vanish a.s K -^ Kc- Accounting for a generic critical 
point at which both the AF and VBS orders vanish is 
in fact highly non-trivial — accoinplished in the theory of 
deconfined quantum-criticality [ij. For the J-K model, 
our QMC results suggest instead that the spin stiffness 
remains non-zero as K —f K~ , and hence that AF order 
is robust in the presence of diverging VBS correlations. 

For K > Kc the VBS states |0)i23 approach the ground 
state exponentially as A^ — > 00. There is thus an N 
above which the VBS quadruplet falls below the lowest 
rotor state |l)o. The rotor states |1,2, . . .)o should also 
develop VBS correlations as i^ ^ i^~, breaking the Z4 
symmetry for K > Kc- Therefore, the rotor states should 
be approached by VBS-like states |1, 2, . . .)i23, which be- 
come degenerate at Kc- The tower of rotor states, along 
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FIG. 4: (Color online) (a) Conjectured behavior of the spin 
stiffness and the spin gap at the AF-VBS transition. (b),(c) 
Possible phase diagrams in an extended parameter space 
where a point with continuous ps can be reached. On the 
solid (blue) curves, the transition is of the type shown in (a), 
with a stiffness jump. The dashed (red) curve in (b) could be 
a line of deconfined quantum-critical points. In (c) there is 
an intervening coexistence phase. 



with the states joining them in quadruplets, thus evolve 
into a series of VBS-like states. This will also apply to 
the fc 7^ towers. What we propose is that these are the 
elementary excitations of the VBS, which hence cannot 
be regarded as "pure VBS" excitations but involve spin 
as well. Clearly there will also be 5* = excitations of the 
VBS — these are the pure VBS excitations, which in our 
scenario would be at higher energy and irrelevant at the 
transition [in contrast to deconfined quantum-criticality, 
where there are gapless 5 = excitations due to an emer- 
gent C/(l) symmetry]. The key aspect of our scenario is 
that all the rotor states change qualitatively (but contin- 
uously) in a manner similar to the ground state, turning 
into gapped VBS/spin excitations. No other low-energy 
states emerge. With no tower of degenerate rotor states 
as iV — > oo there can be no long-range AF order and 
hence the stiffness jumps to zero when the gap opens. 

What is the nature of such spinful gapped VBS ex- 
citations and, specifically, why should they be gapped? 
Regarding a VBS dimer as two spins i,j in a singlet 
(ti-Lj — iiti); the state |l)o can be obtained by exciting 
a dimer into a state (tUj + iitj)- This pair state can 
separate spatially into two spin-i degrees of freedom — 
spinous — between which a string of out-of-phase dimers 
forms |l2l |. The string provides a linear confining po- 
tential for the spinous and hence there should be a dis- 
crete set of gapped excitations. This singlet-triplet pic- 
ture is strictly based on SU{2) symmetry, but even in the 
XY case, where an excited dimer is not a triplet state, 
there should be analogous confining string states. At the 



critical point the spinous will become deconfined, and 
hence our scenario shares some features with deconfined 
quantum-criticality. There are major differences, how- 
ever; the absence of emergent U{1) symmetry and AF 
correlations that are not critical but '^ e""^'', where A 
is the gap and a is a constant. The AF correlation length 
is thus divergent as ii" — *■ Kc from above, but at Kc the 
correlations are long-ranged, not power-law. 

We have not explained why the VBS forms, but we 
have argued that if it does form continuously a simultane- 
ous AF transition of the type we have found numerically 
is not an unlikely consequence. We also cannot predict 
the universality class, because our mechanism does not 
require a particular scaling of the gap at the critical point 
(related to the dynamic exponent). 

We have argued for the existence, in an extended pa- 
rameter space, of an AF-VBS transition in which also the 
spin stiffness vanishes continuously. In Fig. ^b) we con- 
sider the evolution of the transition point when a suitable 
coupling V is added to the J-K model. A line of critical 
points — which would presumably be in the deconfined 
class — could extend from the point at which the stiff- 
ness first becomes continuous. Another possibility is that 
the transition is first-order beyond this point. Fig. El^b) 
shows a scenario where there is coexisting AF and VBS 
order, which requires a different mechanism than the one 
we have discussed. 

In principle we cannot rule out a weakly first-order 
AF-VBS transition 0, [l3 in the J-K model, but our 
results do not require this. Our simulations, in combina- 
tion with general arguments of symmetry-breaking and 
adiabatic evolution of quantum states, instead point to 
a new potential route to quantum phase transitions be- 
yond the Ginzburg-Landau- Wilson framework. To con- 
firm this scenario, it would clearly be important to iden- 
tify a field theory exhibiting this type of transition. 
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